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We present universal relations for a two-dimensional Fermi gas with pairwise contact interactions. 
The derivation of these relations is made possible by obtaining the explicit form of a generalized 
function - selector - in the momentum representation. The selector implements the short-distance 
boundary condition between two fermions in a straightforward manner, and leads to simple deriva- 
tions of the universal relations, in the spirit of Tan's original method for the three-dimensional 
gas. 
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I. INTRODUCTION 

The physics of strongly interacting quantum many- 
body systems has been one of the most active fields 
of research for a number of decades. Strongly corre- 
lated states of matter are ubiquitous in many areas of 
physics ranging from atomic, molecular, and optical 
condensed matter and the study of quantum phase 
transtions [s^, quark-gluon plasma to the physics of 
neutron stars [5]. A prominent example is the observed 
non- Fermi liquid behavior in high-tcmperature supercon- 
ductor materials suggesting a strongly-coupled state be- 
yond the Landau paradigm [gj. Perturbative approaches 
to strongly-coupled quantum systems are untenable and 
even modern numerical techniques can be hard to apply. 
It is therefore of great interest to have analytical insights 
into systems where interactions are strong. 

Research in strongly interacting three-dimensional 
(3D) Fermi gases has experienced an intense growth since 
the recent appreciation of Tan's universal relations Q- 
These important results relate many of the many-body 
properties, such as the adiabatic change of the energy 
when varying the two-body scattering length, the asymp- 
totic momentum distribution, the two-body loss rate [8| 
and the pressure to a single quantity called the contact, 
and were recently verified experimentally • These 

universal relations and quantities involved therein have 
received much recent attention [ll|-|2li. 



Interestingly, Tan relations have also been obtained 
for the one-dimensional (ID) Fermi gas (22j . by means 
of the operator product expansion |23| that had been 
successfully applied in 3D 8]. However, less is known 
about these universal relations for the two-dimensional 
(2D) Fermi gas. Only the energy theorem [23 - l26j . which 
relates the energy of the system to the momentum dis- 
tribution, the large-momentum tail of the latter [2^ [2^ 
and the adiabatic theorem [2^ have been so far derived, 
and only in the homogeneous, or untrapped case. The 
problem has been recently investigated numerically in ref. 

In this paper, we derive the universal relations for the 



zero-range interacting spin-1/2 Fermi gas in 2D, with or 
without an external trapping potential. To do so, we use 
Tan's original method involving a generalized function, 
called Tan's selector. We find the explicit form of the se- 
lector in momentum space, in analogy with its 3D coun- 
terpart recently obtained by one of us [l^ . Our approach 
greatly simplifies the derivation of the adiabatic theo- 
rem with respect to more conventional methods. Once 
the adiabatic theorem is shown, we use it to derive the 
generalized virial theorem, the pressure relation and the 
two-body loss rate. 



II. SYSTEM HAMILTONIAN. 

The zero-range interacting many-body Hamiltonian for 
N = N-f + Ni spin-1/2 fermions has the form 
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where m is the particle's mass, r^j = | 
single-particle external potential, and V is the 2D regu- 
larized pseudopotential [2^ 
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which is a particular member of the Olshanii- 
Pricoupenko A- family of pseudopotentials [s^l • Note that 
pseudopotential ([U does not depend on the 2D scat- 
tering length a, which only enters the problem through 
the boundary condition at short-interparticle distances 
i^i^ij) ^ \og{rij / a). The universal singlet binding en- 
ergy in the homogeneous case {W = 0) can be seen to be 
[31 
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where 7 = 0.577215665... is Euler's constant. The bind- 
ing energy is the most appropriate fitting parameter of 
the theory in momentum space, as we see in the following 
paragraph. 
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III. TAN'S SELECTOR. 

To derive the universal relations for the Fermi gas in 
2D, we will make use of Tan's original method 0], par- 
ticularized to two dimensions. In ref. 0], Tan de- 
fined two selectors, fj and Z, the last of which is un- 
necessary in 2D, and involves the logarithm of a dimen- 
sioned quantity. The ry-selector satisfies Jd^rr)(r) = 1, 
/ d^rfi{r) log(r/a) ~ 0, while fi{r) = for r 7^ 0. 

We here redefine the momentum representation r](k) 
of ry(r) as follows 



r]{k) = 1 (fc < 00), 



fPk'^/m + EB 



0, 
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where Eb is the two-body universal binding energy, Eq. 

It is easy to see that Eq. ([5]) is equivalent to Eq. 
(16b) of ref. 2^, but our condition is more natural and 
represents a clear way of writing the two-body bound 
state integral equation. Proceeding in analogy to the 3D 
case 128il, we see that 



where R' = (r2, . . . ,rAr^, S2, . . . ,sn^), r = ri - Si and 
ro = (ri -I- Si)/2, and K is defined as 

if(R',ro,t)= f d2r0*(R',ro,r)0(R',ro-l-t/2,r + t), 

Jr<e 

(10) 

where 0(R',ro,r) stands for (/>(R, S), with ri and si re- 
placed by ro-l-r/2 and ro — r/2, respectively. In Eq. (|10p . 
e > is a small real number, however much larger than 
t = |t|. The goal is to show that the integral in Eq. © 
vanishes since, if it does, relation ([7]) holds [7.]. First, we 
expand 4> as 



(R', ro, r) - A(R', ro)[log(r/a) + 0(r)], 



(11) 



and 



A{R', ro + 1/2) = A{R', ro) -t- Vr^ A(R', ro) • t/2 + 0{f). 

(12) 

We obtain K = Kq + Ki , with 
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The above explicit form for r](k) will prove crucial for the 
derivation of the two-dimensional adiabatic theorem. 



IV. ENERGY THEOREM. 

In ref. [IJ] Tan stated, without derivation, that the 
energy of the 2D Fermi gas can be written as 
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where riko- — {c^ „c\^.a): with Ck.o- the spin-cr annihilation 
operator at momentum k, is the momentum distribution 
and W = the total external potential. 

To prove Eq. ^ it is sufficient to consider a pure, 
normalized state |0) with No- fermions of spin a =tii 
and total number of fermions N — + Ni 



dRdS0(R, s) n ^t('"^) n '/^i^^^-) 
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where R (S) is shorthand for ri,...,rAr^ (si, . . . , sat^) 
and the integrals are done over the whole 2A^-dimensional 
space. Above, ipti^) the spin-a creation operator at 
position r, and |0) is the vacuum. 

Tan's derivation Q for the 3D gas can be followed in 
parallel for the 2D case until the following expression is 
encountered 



dR'd^ro j dHfi{t)\7lK {R',ro,t), 



(9) 



K, « 7rA*(R',ro)-(log(Va) - 1) A(R', ro) • t. 

(14) 

Inserting K into (|9]), we see that the integral vanishes, as 
we wanted to show. 

We can express the energy relation in Eq. ([7]) more 
explicitly by using the form of the ry-selector in Eq. ([6]) , 

as 
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where $7 is the area and C — linik->.oo k^ny^.a is the 2D 
contact density. The above relation, Eq. (fT5|) . was ex- 
pressed in a somewhat different manner by Combescot et 
al. in [25i] and by Werner and Castin in |26|. The tail of 
the momentum distribution was found in |25l [26j to be 
oc 1/fc^ and therefore the contact C is a finite quantity. 



V. ADIABATIC THEOREM. 

The adiabatic theorem states the relation between the 
contact and the change in energy as the scattering length 
is slightly varied. This is the central result among the 
universal relations for the Fermi gas, since it is needed to 
derive the virial theorem, the pressure relation and the 
inelastic two-body loss rate. 

The fact that the 2D pseudopotential, Eq. ([2|), does 
not depend explicitly on the scattering length adds a 
complication to the proof of the adiabatic theorem that 
was not present in 3D and ID. Indeed, the Hellmann- 
Feynmann theorem, which was employed to derive this 
important result in 3D by Braaten and Platter [8] and in 
ID by Punk and Zwerger [l^], cannot be applied directly 
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to Eq. ([2]) nor to any member of the so-called A-family 
of DseudoDotential s f29l Isoj . In addition, Tan's approach 
to the 2D problem [23| lacked a simple dependence of the 
ry-selector on the scattering length. These two facts have 
prevented a simple derivation of the adiabatic theorem 
using regularized pseudopotentials, which has been first 
shown by Werner and Castin [26^] who obtained it using 
Bethe-Peierls boundary condition. Below, we present a 
very simple proof of the adiabatic theorem by using the 
explicit form, Eq. of the 77-selector. 

We begin by defining the following operator 



(16) 



where ry(fc) corresponds to scattering length a. It is then 
easy to see that 



E{a') - E{a) 



{cp{a)\r^aT-r^a,T\^{a')) 
{4>{a)\4>{a')) 



(17) 



where E{a) and 4>{a) are an energy eigenvalue and as- 
sociated eigenstate of Hamiltonian ^ corresponding to 
scattering length a. Using Eq. ([6]) wc find 

E{a') - E{a) _ ^ h'^P6{l/k) (0(a)| 4 ,,Ck,a |0(a')) 



log(a7a) 



E 



2mk 



(0(a)|0(a')) 



(18) 

If we now take the limit a' — a in the above equation, 
we obtain the desired adiabatic theorem 
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An immediate application of the adiabatic theorem is 
the calculation of the contact for a weakly-coupled 2D 
Fermi gas. The energy density for the spin-balanced case 
is given by [s^l 



E 
N 



Am 
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1 



log(fcFa) 
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where N is the number of particles and kp is the Fermi 
momentum. The contact in this limit is therefore given 
by 



C = 



npkj, 



2\og^ {kpa)' 
where p = N/ft is the particle density. 



(21) 



VI. GENERALIZED VIRIAL THEOREM. 

We assume now that the system is under an external 
potential of the form W{r) oc , which has not been 
considered before for the 2D case. In this case, we will 
show that the virial theorem reads 



E : 



Attiji 
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To see this, we follow a technique used for the 3D gas 
0] . We begin with a state 4> corresponding to scattering 
length a. We slightly change the scattering length to 
a' = {l + e)a, with e small, associated with state 0' which 
has energy E' . Expanding E' in powers of e, and using 
the adiabatic theorem, Eq. (fT9|) . wc find 



E' ^ E 
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We now define a rescaled wave function (/)"(R, S) = (1 + 
e)^(/)'((l + e)R, (1 + e)S). Expanding 0" as in Eq. ([TT]), 
we find that 4>" is a state at scattering length a" = a'/(l+ 
e) = a. Using the scaled wave function in the expectation 
value of the energy, we find its energy E" to be 

E" - (1 + efEi, + (1 + ,)-P{(p'\ W 10'), (24) 

where E[^ — E' — {4>'\W Expanding E" in powers 
of e, we obtain E" - E' = 2eE[^ - /3e((/)'| W \4>') + 0{e^). 
From the quadratic convergence properties of variational 
energies [34!| we have E" — E — O(e^), and therefore, 
using Eq. ([23t twice, we obtain 



2'Kra 



0{e'), (25) 



which, after taking the limit e — >■ 0, proves the virial 
theorem (l22l). 



VII. PRESSURE RELATION. 

For a homogeneous system {W = 0), there is a relation 
between pressure, energy and contact, which has been 
derived in 3D 0,[i] and ID '^]. In 2D, it is given by 



P = 



E h^C 



(26) 



To show Eq. ((26)) . we follow most of the steps taken for 
the proof of the virial theorem ([22]) , and we arrive at the 
expression 



pn = hm 



E"-E 



1 - (1 + e)- 
27rm 



2E 



lim 



>ol-(l + e) 



(27) 



Using L'hopital's rule in the above equation, we obtain 
the desired pressure relation ([^5)) . 



VIII. INELASTIC TWO-BODY LOSS RATE. 

Two-particle loss rates can be calculated |8| by adding 
a small imaginary part a/ (|a/| ^ \a\) to the scatter- 
ing length, so that we replace a — > a + ia/. The re- 
sulting complex energy is then expanded as E{a + iaj) = 
E{a) + iaidE /da{a) + 0{{ai)'^), and the inelastic loss rate 
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r is identified as —T/2 = ajdE/da. Using the adiabatic 
theorem p^ . we obtain 

r = -a,^ + 0((az)2). (28) 



IX. CONCLUSIONS. 

We have obtained the universal Tan relations for the 
two-dimensional spin- 1/2 Fermi gas with zero-range in- 
terparticle interactions. These results can be tested with 
current experimental techniques such as Bragg [13] or RF 
spectroscopy d, IH, Isst, already available for the two- 
dimensional Fermi gas Our methodology is close 
to the original approach for the three-dimensional prob- 
lem The power of our approach relies on obtain- 
ing the explicit form, in momentum representation, of 



the so-called Tan's selector, which expresses the contact 
conditions between two particles with a given scattering 
length in closed form. This allowed us to overcome the 
fact that the Hamiltonian of the system does not depend 
on the two-body scattering length, and led us to a more 
straightforward derivation than with conventional meth- 
ods. Our method is directly applicable also to univer- 
sal three-body physics in two dimensions (37l - [39j where 
additional three-body contact parameters were recently 
proposed in the three dimensional case [13, |4l| . This is 
an interesting direction for future study. 
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